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Abstract 

Hard thermal loops play a central role in the theory of long wavelength excita- 
tions of a quark-gluon plasma. We show in this paper how our recent derivation of 
their generating functional from the Dyson-Schwinger equations sheds light on their 
semi-classical nature and their remarkable gauge structure. In particular, we show 
that our kinetic equations can be written in terms of auxiliary gauge fields with 
zero curvature. Remarkably, the latter property determines entirely the dynamics 
in the kinematical regime of hard thermal loops. This explains in particular why 
the generating functional could have been obtained by Taylor and Wong from gauge 
invariance considerations. It also clarifies the role of Chern-Simons theory in this 
context. 



*CNRS 

^Laboratoire de la Direction des Sciences de la Matiere du Commissariat a l'Energie Atomiquc 



1 Introduction 



"Hard thermal loops" (HTL) are the dominant contributions at high temperature to n- 
point amplitudes with soft external momenta (P ~ gT, with T the temperature and 
g the gauge coupling, assumed to be small). Their role in perturbative calculations, 
and their remarkable properties, where first established via a systematic diagrammatic 
analysis^, |], ||]. One of their most noticeable property is their gauge symmetry; specif- 
ically, hard thermal loops obey Ward identities similar to those of the tree amplitudes, 
suggesting that their generating functional should be a gauge invariant functional of the 
soft fields. The gauge invariance has been used by Taylor and Wong to construct a gen- 
erating functional for the HTL's in closed form[Q. Later, this effective action has been 
rewritten in a manifestly gauge invariant form by Braaten and Pisarski[Q and also by 
Frenkel and Taylor @. It was also realized that the generating functional of HTL's could 
be understood as the eikonal of a Chern-Simons gauge theory, raising the hope that the 
use of the known properties of the latter theory could lead to a better understanding of 
gauge theories at finite temperature 0. 

In a series of recent papers, we have shown that the hard thermal loops are intimately 
related to the long wavelength, collective, excitations of the hot gauge plasmas, much 
alike the non trivial dielectric properties of the electromagnetic plamas are related to the 
collective oscillations of their charge constituents. To establish this correspondence, we 
have developed a theory of the collective excitations, which is gauge covariant and takes 
consistently into account all the leading terms in an expansion of the equations of motion 
in power of the coupling strength || ^ [TO]. 

A central quantity in our approach is the "induced current" which represents the 
response of the plasma particles to applied background gauge fields. This is obtained, 
after solving nonlinear equations which may be given the form of kinetic equations, as a 
non local functional of the gauge fields. The nonlinearity of these equations is entirely 
due to gauge rotations of the current (a vector in color space) induced by covariant 
derivatives. The non locality of the induced current arises from the need, in calculating 
the response functions, to integrate the history of the hard particles along their straight 
line trajectories. These simple features, whose physical origin is tranparent, explain most 
of the properties of the HTL's, properties which remained somewhat surprising in the 
diagrammatic approach. 

In this paper we wish to pursue our analysis, and make closer contact with the 
formal works mentioned above. In particular we shall show that the gauge structure of 
the HTL's follows from a simple property of our kinetic equations, which is ultimately 
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related to a property of the Lorentz equation describing the motion of a classical particle 
in an electromagnetic field, or its generalisation to the non abelian case. More precisely, 
we shall show that, in the same way that the Lorentz equation can be rewritten in terms 
of the gauge dependent canonical momentum, our kinetic equations can be also written 
in terms of auxiliary gauge potentials with zero field strength. The latter property follows 
from the fact that, within the present approximation scheme, the plasma particles do not 
change direction in the course of their interaction with the gauge field. Quite remarkably, 
in this kinematical regime, the condition of zero field strength turns out to be sufficient 
to determine the dynamics. As we shall see, this provides an explanation for the fact that 
it has been possible to construct the generating functional for the hard thermal loops on 
the basis of gauge invariance only, with no dynamical input beyond the form of the gluon 
self energy M, 

The plan of the paper is as follows. In the next section, we present a short summary 
of our approach, and recall the form of our basic kinetic equations. In section 3 we discuss 
some of their semi-classical features, and show that they can be given a simple gauge 
structure. Then in section 4 we show explicitely how the gauge invariance condition used 
by Taylor and Wong to derive the generating functional of HTL follows from our kinetic 
equations. We also clarify the role of Chern-Simons theory whose relevance in the present 
context has been much advocated recently. Finally, we verify that the expression for the 
energy density obtained by Nair using Chern-Simons theory coincides with the one that 
we have obtained by integrating the mean field equations. The last section summarizes 
our conclusions. 



2 Semiclassical theory of collective excitations 

Our semi-classical theory of the collective excitations is based on a perturbative anal- 
ysis of the QCD Dyson- Schwinger equations. It is obtained by consistently preserving 
the leading terms in an expansion in powers of the gauge coupling g. In doing so, one 
encounters three types of approximations, which, in most many-body systems, are inde- 
pendent approximations, but here are controlled by the same small parameter, i.e. g. 
These are the weak coupling approximation (g <C 1), the long wavelength approximation 
(A ~ 1/gT ^> 1/T), and the small amplitude approximation (the gauge field strength 
tensor is limited by F ^ gT 2 , or, equivalently, the gauge potentials satisfy A&T). 

The result of our analysis is a set of coupled equations. The first equation, the gen- 
eralization of the Maxwell equation in a polarizable medium, relates the gauge mean field 
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A a ^ to the induced current j't (see (a) below). The other equations allow the calculation 
of the induced current in terms of the average gauge fields. They can be given the form of 
simple kinetic equations, analogous to Vlasov's equation of ordinary plasma physics (see 
(c) below). In our approach, the hard thermal loops appear as functional derivatives of 
the induced current with respect to the gauge fields || |9|, . 

(a) Mean field equation 

The equation of motion for A ^ takes the form 



[D»,F v »{x)] a = (2.1) 

where = <9 M + igA^x), (A^ = Ap a ), F^ = [D^D v \j(ig) is the field strength ten- 
sor, and j)j,{x) is the induced current, which transforms as a color vector in the adjoint 
representation. (The generators of the gauge group in different representations are taken 
to be Hermitian and traceless. They are denoted by t a and T a , respectively, for the fun- 
damental and the adjoint representation, and are normalized so that Tr(t a t 6 ) = (l/2)5 ab 
and Ti(T a T h ) = N5 ab .) 

(b) Structure of the induced current 

The induced current expresses the response of the plasma particles to the long 
wavelength color fields A a ^. It is convenient to separate the contributions of the fermions 
from those of the transverse gluons, by writing j" = j f a M +Jw For the fermions (with Nf 
flavors) we have 

Jf ;(x) = gN f j-0-v IM (8nl(k,x)-5n^,x)), (2.2) 

where = (1, v) and v = k/k is the velocity of the hard particle (k = |k|, and |v| = 1). 
The quantities 5n^_(k, x) are the components of a color vector in the adjoint representation. 
They may be combined into a density matrix, 5n± = 5n± t a , ultimately related to the 
fermion propagator in the presence of the gauge fields. There is no color singlet component 
in 5n± because we are considering only excitations which carry gluon quantum numbers. 

Similarly, the current carried by the transverse gluons can be written in terms of a 
density matrix 5N(k, x) = 5N a (k, x) T a : 



29 1 (^^TrT^iV(k,x) =2^iV J ^v,5N a (k,x). (2.3) 



(c) Kinetic equations for Sn± and SN 
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The color density matrices 8n± and 5N are determined by the following kinetic 
equations || |T0 



[v ■ D x , 8n±(k,x)} = T9 v • E(x) (2.4) 

[v ■ D x , 8N(k, x)] = -gv E(x) (2.5) 

where E l a = is the chromoelectric field. Furthermore, Nf. = l/(exp(/?fc) — 1) and 
n fc = l/(exp(pk) + 1) are, respectively, the equilibrium boson and fermion occupation 
factors. Note that, in deriving these equations, we have used the isotropy (in k space) of 
the equilibrium distribution functions. For more general, anisotropic, equilibrium distri- 
butions n(k) and iV(k), we should replace 

in the r.h.s. of eq. fl2.4|) , and similarly for the second equation. 

In the abelian case, eq.(2.4) reduces to the well known linearized Vlasov equation 
describing the evolution of single particle distribution functions. It is perhaps useful here 
to emphasize that, although the underlying dynamics is that of the Lorentz equation 
describing the motion of a classical charged particle in an electromagnetic field (see next 
section), however 5n(k,x) is a Wigner function whose only slow variations in x have a 
classical character (in particular, the plasma particles have typical wavelength comparable 
to their mean separation and remain therefore quantum particles). In the non abelian case, 
the distribution functions are matrices in color space: the color degree of freedom remains 
essentially quantum. (In this respect we differ from Heinz who has given a derivation of 
equations similar to ours starting from classical kinetic theory, and treating color also as 
a classical degree of freedom [Til].) 



For the subsequent analysis of these equations, it is convenient to express 8n± and 
8N in terms of new functions, W^{x; v) = W£(x\ v)t a , defined as the solutions to 

[v ■ D x , W»{x-v)\ = F» u (x)v„. (2.6) 

The quantities W£(x; v) are the components of a color vector in the adjoint representation, 
and satisfy 

v^(x;v) = 0, (2.7) 
so that = v % W % a . From the equations above, it is easily seen that 

5n a ± (k, x) = T 9 W°(x; v) ^, 5N a (k, x) = -gW° a (x; v) ^. (2.8) 
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Note that for an anisotropic initial distribution, the fluctuations 5n± and 5N a depend on 
all the three independent components of W^(x; v), that is, e.g., 

On 

6n%{k,x) = TgW l a {x-v)— v 

and similarly for SN a . 

(d) Solving the kinetic equations 



The equation ( |2.6|) for W£ is easily solved once boundary conditions are given. 
We assume here retarded conditions, such that the fields vanish as io — ■> — oo. The 
corresponding solution of eq. (|2.6|) reads 



W(x;v) = / drU(x,x - vr) F^(x - vr)v u U(x - vr,x), (2.9) 
Jo 

where U(x, y) is the parallel transporter along the straight line 7 joining x and y, 

U(x,y) = Pexp{-ig f dz^A^z)}, (2.10) 
with P denoting the path-ordering operator. Thus, 

U(x,x — vt) = P exp j— J dsv ■ A(x — v(r — s))\ . (2-H) 



In order to verify that Q2.9Q is the correct solution to eq. ( p.6|) we may use the following 
formula for the line-derivative of the parallel transporter 

(v-d x )U(x,y) = —igv ■ A(x)U(x,x — vt). (2-12) 

y=x—vr 

We shall use later the quantity 

V(x:v) = lim U(x,x—vr), (2.13) 
which is well defined since A^x — vt) — > as r — > 00. We have 

{y ■ d) V(x; v) = -igv- A(x) V(x; v), (2.14) 
which leads to the following representation for v ■ A(x): 

v ■ A(x) = -((v d)V) V" 1 . (2.15) 
9 v ' 

The matrix V(x; v) generates the gauge transformation from the gauge v ■ A = to the 
actual gauge. 

(e) Induced current in closed form 
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Once the solution of the kinetic equation is known, one can easily calculate the 
induced current in closed form. By inserting eqs. (|2.8|) in the expressions ( |2.2|) and (|2.3|) , 
we obtain 

dVL 



3w 



Alt 



v"W^(x;v). 



(2.16) 



after the integration over k — |k| has been performed. Here u v is the plasma frequency, 
ui = (2N + Nf)g 2 T 2 /18. The integral / dfl runs over all the directions of the unit vector 



v. According to eqs. ( [2.9|) and (|2.16| ), the retarded current reads 

dVL 



f{x) 



3ul — v» dr U(x, x — vr) v • E(x — vr) U(x — vr, x). (2.17) 
i w io 



This expression can be viewed as a generating functional of the hard thermal loops with 
retarded boundary conditions. It was first presented in Ref. ||, and has been rederived 
since using other methods [p"3| (see also below). The current ( |2.16|) is covariantly conserved, 



D»,f{x) 



0. 



(2.18) 



as necessary for the consistency of the mean field equation ( |2.1|) . The "abelian-like" Ward 
identities relating the HTL's, which were first identified in the diagrammatic approach of 
Refs. 0, 131 direct consequence of eq. 



3 Gauge structure 

We turn now to a physical interpretation of some of the quantities that we have introduced, 
and to an analysis of the simple gauge structure which emerges from our equations. 

First, it is useful to recall various forms of the equations of motion for a classical 
particle of mass m and charge e, moving in an electromagnetic background field. From 



the Lagrange function[12 



L = — mVl — v 2 — ev ■ A, (3-1) 
one easily obtains (t> M = (1, v)) 

t = - e "" v - 4 - (3 - 2 » 



The canonical momentum 



dL mv A . 

p = — = -== + eA = k + eA 3.3 
<9v VI - v 2 
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is the sum of the kinetic momentum k, related to the velocity of the particle, and e times 
the gauge field A. The energy of the particle may be given a similar form, 



TO 

p° = p • v - L = - + eA° = k° + eA°. (3.4) 



It obeys the equation 



dv° 

= evPdoAv (3.5) 



so that, for the 4- momentum we have: 

% = ev v d»A v . (3.6) 
at 

This equation is not manifestly gauge covariant. However, the equation for the kinetic 
4- momentum k^ = p 11 — eA 11 , 

dk^ 



eF^(x)v l/ , (3.7) 



dt 

is explicitly independent of the choice of the gauge. We shall see that, to the two ways of 



writing the Lorentz equation, namely eqs. (|3.6|) and ( |3.7|) , correspond two ways of writing 
our basic kinetic equation ( j2.6|) for W^. 

Consider first a QED plasma, where eq. Q2.6|) reduces to 



(vd x )W^(x;v) = F^(x)v u . (3.8) 

This is very much similar to eq. ( p.7| ). To make the analogy closer, we remark that the 
line derivative in the l.h.s. of eq. (|3.8| ), i.e. v ■ d x — d t + v • V, may be interpreted as the 
total time derivative along the trajectory x(t) = x + vt of a fictitious particle. That is, 
eq. ( |3.8| ) may be rewritten as 

±W»(t,x{t);v) = F^(t,x(t))v u , (3.9) 

showing that edW^(x\v) is the kinetic 4-momentum acquired during the time dt by a 
charged particle due to its interaction with the electromagnetic field. Note, however, that 
in calculating according to eq. ( |2.9| ) we have assumed the velocity v to be constant, 
while in the Lorentz equation (|3.7|) , k^ and v are related quantities (fc M = kPv^). This is 
in line with our assumption that the typical momenta of the thermal particles are "hard" , 
k ~ T, while those of the background fields are "soft", P ~ gT k. Thus the hard 
particles follow straight line trajectories at the speed of the light, and these trajectories 
are not altered, to leading order, by the interaction with the background gauge field. The 
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condition ( f2.7| ) reflects then simply the fact that the energy transferred by the field, eW°, 
coincides with the mechanical work done by the Lorentz force, ev • Ak = ev l W l . With 
this interpretation, the formulae (j2.8[) for 5n± become transparent: 



* !± (k,x, ( ) = ^f^W'^-^ 

« n(k- Ak±(x,*)) - n(k), (3.10) 

that is, the particles found at time i at the point x with momentum k, are those which, in 
the initial distribution, have momentum k — Ak(x, t), where Ak(x, t) is the momentum 
acquired by the particles whose trajectories go through x at time t. 

In the non-abelian case, the fluctuations Sn± and SN are matrices in color space, 
that is, = W£t a . The color vector of components W£ precesses in the background 
gauge field. This precession is induced by the covariant derivative in eq. (|2.6| ). Viewing 
this precession as an additional source of time-dependence for the color vector one 
can write 

j t W»{t,*{t)-v) = [(d t + v V)5 ac - gf abc (v-A b )]w?, (3.11) 

so that eq. ( |2.6| ) may be given a form similar to eq. ( |3.9| ). 

We shall see now that it is possible to rewrite our kinetic equation in a form similar 
to eq. ( p.6|) for the canonical momentum. This follows by noticing that 



Thus, if one sets 



one gets 



Ft*" v u = d»{v -A) - v ■ D, A» . (3.12) 



a"{x; v) = A"{x) + W{x; v), (3.13) 



[v-D, a"] = d^(v ■ A), (3.14) 
showing that is a functional of v ■ A only. Explicitly, 

a^{x;v) = I drU(x,x — vr) d^[y ■ A(x — vr)^j U(x — vr,x), (3.15) 



where the retarded boundary conditions have been taken into account. In the abelian 
case, eq. ( |3.14j ) can be rewritten in a form analogous to eq. ( |3.6| ) for p^: 



Y a»(t, x(t); v) = {d t + v • V) = d"{v-A), (3.16) 

LLC 



S 



showing that go? may be indeed understood as the change in the canonical momentum 



+ gA^ 1 of a fictitious particle following the trajectory x(£) = x + vt. 



The two writings of our kinetic equations, namely eq. ( |2.6| ) and eq. ( p. 14 ), have their 
respective advantages. Eq. ( |2.6|) is manifestly gauge covariant, while eq. ( |3.14|) is not. On 
the other hand, the new fields a^ 1 transform as gauge potentials, which has interesting 
consequences, as we shall see now. Under a gauge transformation induced by the operator 
h(x) = exp(—ig9 a (x)t a ), = W£t a transforms as W fl (x) — > h(x) W^(x) h 



-l. 



x 



while 



A, 



A-t a 



hA^h' 



- (ij 'g)hd tl h~ 1 . Thus a^(x;v) transforms indeed as a gauge 
potential. Let T> M = <9 M + ig be the covariant derivative constructed from a M , and 
fnv{x',v) = [D^,V u ]/(ig) the corresponding field strength tensor. Then it may be shown 
that 



= du,a u - d v a„ + ig [a^, a v ] 



(3.17) 



To verify ( |3.17| ), one could directly derive the equation satisfied by f^, starting with 
eq. ( |3.14| ). This is the way followed in Ref. [jl? , where eq. ( 3.17 ) was first obtained (see 
eq. (3.24) in Ref. |17]])- An alternative derivation is to solve eq. ( 3.14j) for a M , with v ■ A 
given by eq. ( |2.15| ). It is then easily seen that that 



a^(x\ v) 
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(3.18) 



from which eq. ( |3.17| ) follows. 

In summary, our basic dynamical equation can be written either in terms of 
(see eq. (|2.6| )), or in terms of a M (as in eq. Q3.14 )). In the abelian case, this freedom 
corresponds to the choice of writing the Lorentz equation in terms of the kinetic or of the 
canonical momentum, respectively. The fact that the field a M turns out to be a pure gauge 
(see eqs. (|3.17|) and (]3.18|) ) results ultimately from the assumption that the trajectories 
of the plasma particles are not deviated by the interaction with the gauge fields. In 
particular, a M can be set equal to zero by an appropriate choice of the gauge, specifically, 
by choosing the gauge v- A — (see eq. ( |3.15|) ). Then, W^ 1 = —A^ 1 , as can be seen directly 
on eq. fl2.6|). Note, however, that for practical purposes, such a gauge is of no use, since 
the calculation of the induced current — the quantity one is primarily interested in - 
involves an integration over all directions of v (see eq. ( 2.16[) ). 

Quite remarkably, the gauge structure that we have exhibited completely determines 
the underlying dynamics. This explains why most of the dynamical information could have 
been reconstructed from arguments based on gauge invariance only[H|. To see this, we 



note that the dynamical equation ( |3.14j ) for a M = A^ + follows from eq. ( |3.17| ), once 
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the condition (|2.7| ) is taken into account. Indeed, by multiplying eq. ( |3.17| ) by v u , we 
obtain 

{v ■ d)a^ — v y (d^a u ) + ig[v ■ a, a M ] = 0. (3.19) 

Using the fact that v M is a constant vector, together with the property v ■ a = v ■ A (which 
follows from eqs. (|2.7|) and (|3.13| )), we can write this equation in the form (|3.14j) . To 



the extent that one is interested only in the calculation of the induced current — which, 
according to eq. Q2.16| ), depends on W° = a — A — only the component \i = of 
eq. ( |3.14j) is needed. Explicitly, 



(v ■ d)a - d (v ■ A) + ig [v ■ A, a ) = 0. (3.20) 

This is, within minor notational changes, the equation which has been obtained by Taylor 
and Wong as the condition of gauge invariance of the effective action |4j]. This is discussed 
in detail in the next section. 



4 Relation to previous works 

As we have seen, in the computation of the induced current, the important gauge variables 
are A and v ■ A = A — v • A. (Recall that a is a functional of v ■ A only, as implied by 
eq. ( p.20|) . ) Alternatively, and in order to facilitate the comparaison with previous works, 
it is convenient to use the variables A + = v ■ A and A_ = v' ■ A, where v M = (1, v) and 
v ,fl = (1, — v). We also set d + = v ■ d and <9_ = v ' ■ d. In these notations, eq. ( |3.20|) reads 

d + a - d A + + ig[A+, a Q ] = 0. (4.1) 

We show now that this is precisely the condition for the gauge invariance of the generating 
functional of the HTL's, denoted by S'fA^]. To do this, let us recall that the induced 
current is, formally, the first derivative of the effective action S, 

dS jS(x) = 3ut[^-v>>W°{ar,v). (4.2) 



5A*(x) Jay J P J 4vr 



Here, "formally" means that in writing eq. ( |4.2| ) we are renouncing to the retarded con- 
ditions for A^ and j M . Indeed, eq. ( |4.2j ) implies 



5Al(y) 5A%(xY 1 ' ) 
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which can only hold on the space of fields in which the operator v -D is invertible|4], |I0| . In 
this functional space, which we denote by 71, we can integrate ( |4.2|) to get S[A\. Setting 
W° = a — A , and separating the contributions of a and A , we write eq. (|4.2|) as 



5S 



6A*(x) 



(4.4) 



where we have indicated explicitly that a depends only on A + = v ■ A. The effective 
action is therefore of the form: 



S = ^ljd^- l -A° a Al + /^W[A + ]| 



(4.5) 



where, by construction, the functional W satisfies 

6W 



a° a (x;v). (4.6) 



SA%(x) 

Consider now an infinitesimal gauge transformation 

Ap^Ap + SAp, 5A^ = [D^6}. (4.7) 
Under such a transformation, S ^ S + 5S, with 

5S = 3 ul J d 4 x j^e a {d A a + - d + a a + gf abc A b + a c } . (4.8) 



Since Oq satisfies eq. (|4.1|) , it follows that 5S = 0. 

The form of the effective action given by eq. ( |4.5| ) was first inferred by Taylor and 
Wong from an analysis of the HTL'sQ. By assuming S to be gauge invariant, they 
obtained eq. (|4.1| ) as an equation for 5W/5A+. And, in fact, our arguments above have 
been intentionally developed in a way which makes the relation to the analysis of Taylor 
and Wong obvious. But, of course, we could directly infer the gauge invariance of S from 
the conservation of the induced current, through a standard, Noether-type, argument. 
That is, by using eqs. (|4.2|) and ( f4.7| ), we can write 

6S= f d 4 xj» 5A fl = - [ d A x 9 a [£>„, A , (4.9) 



where the second equality follows after an integration by parts (the surface term has been 
assumed to vanish). Because of the conservation law for j^, eq. Q2.18D , the last integral 



vanishes, and so does SS. The relation between the condition for the gauge invariance of 
the effective action and the dynamical equation for has also been pointed out in Ref. 

ra. 
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The arguments above suggest a different way to compute the induced current, which 



has been, in fact, followed by Jackiw and Nair in Ref. [13]. Starting with the effective 
action S of Taylor and Wong, eq. ( |4.5| ), they have constructed j% as the first derivative of 
S, i.e. eq. ( |4.4| ), with a°[v4+] = 5W/5A+ determined from the gauge-invariance condition 
on S, eq. (|4.1| ), as in the original work of Taylor and Wong. In order to recognize that 
the expression for j£ obtained in Jl3|] is the same as ours, we follow Ref. and replace 
eq. ( |4.1| ) with an equivalent equation for the new function a_ = 2a® — A + : 



d+ a- — d- A< 



19 [A, 



0. 



Then, the current 



jHx) 



becomes 



Air 



v" (oq - A ) 



3 2 
2^ 



dVL 
Air 



(a. 



A. 



(4.10) 



(4.11) 



We recognize in eqs. ( |4.1U| ) and ( [4.1 1|) the expression of j M proposed in Ref. [|T^]. 

On the other hand, it can be easily verified using eqs. (|2.7j ) and (|3.13|) that a_ just 
defined is also equal to v' ■ a. This is not surprising, since eq. ([4.1 0|) also follows from 



eq. ( p.!4| ) for a M after multiplication with v' . Let us further define a H 



v ■ a. 



Because 



of eqs. fl2.7|) and (|3.13|) , we have a + = A + (i.e., v ■ a 
eq. (|3.19|) ), so that eq. ([4 . 1 0|) may be finally written as 

d + a_ — 9_ a + + ig[a + , a_] - 



A, as already noticed after 



(4.12) 



which coincides with the projection of the zero-curvature condition ( p. 17] ) onto the hyper- 
plane defined by the 4- vectors v M and v that is, eq. ( 4.12| ) is the same as f 



0. 



If we recall that a + = a — v ■ a and a_ = ao + v ■ a, we see that the gauge structure 
defined by eq. ( [4.12j ) involves the temporal (ao) and the longitudinal (v ■ a) components 
of a M , but not its two transverse components. This equation represents the starting point 
of the relation to the Chern-Simons theory explored in Refs. @. 

The main difference between the two derivations of j M outlined above concerns 
the method used for solving eq. ( |4.10| ) (or, equivalently, eq. ( |2.6| ) for W°). We get the 
solution directly in Minkovski space, where retarded conditions are easily implemented 
(see eq. ( |2.9| ) above). The authors of Ref. [0 solve eq. ( |4.10| ) in Euclidian space-time, 
using technics borrowed from the study of the Chern-Simons theory. They obtain thus a_ 
as a power series in A + |7j. In order to get the retarded current, an analytic continuation 



has to be performed, to ensure the correct ie prescript ion [|T^ . 

The functions W 11 have been found useful in the construction of the energy-momentum 



tensor for the gauge fields [0, [17]], and also in the construction of a manifestly gauge invari- 
ant expression for the effective action. Such an expression was proposed first by Braaten 
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and Pisarski0. Their expression reads (with (DF) a = [D,F] a 

3 

/ 

4:71 



s 



2^ 



J d 4 y Tr 



F tJ ,x(x)(x\ 



v-D) 



(4.13) 



and has been rederived since by Frenkel and Taylor [11] (who inferred S from a study of 
the forward scattering amplitudes for quarks and gluons on "soft" gauge fields), and by 
us (in Ref. fl(|], we have shown explicitly that the functional derivative of fl4.13|) gives 
the correct induced current). The functional ( [4.13 ) is well defined only on 71, where it 
coincides with|fT4}| 



~2 U * 



d x 



dVL 

4lX 



TrW< l (x;v)W tl (x;v) 



By using eqs. 



13]) and (|3.18| ), we can write 

% 



W»(x;v) = -(D»v)v-\ 
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Thus, an alternative expression for S is 



S 



3 w 
2^ 



Jd*xJ§Tr(B»v)^v)\ 



(4.14) 



(4.15) 



(4.16) 



which also appears in Ref. [15 



Let us finally consider the energy £ of a gauge field configuration in the plasma. 
Different approaches has been proposed recently to compute this quant ity[[TJ], [15|, 16, 17], 
and the various results obtained are not related in an obvious way. In Ref. fl7| , we have 
calculated £ by integrating the field equations of motion (|2.1|). One of our results may be 
written in the form £ = £ YM + £ ind , where £ym is the standard Yang-Mills contribution, 



£ YM = J d 3 x i(E a (x) • E a (x) + B a (x) ■ B a (x] 
{Bi(x) = -(l/2) e y fc F a ^(x)), while 



£ind — 3co> 



d 3 x I — TyW°(x;v)W°(x;v) 
4tt 



(4.17) 



(4.18) 



may be interpreted as the polarization energy of the plasma, that is, the energy transferred 
by the gauge fields to the plasma constituents JI7J. According to these expressions, the 
field energy is obviously positive. A similar conclusion has been reached by Nair|l6|, who 
derived a Hamiltonian for the soft gauge fields by exploiting the analogy with the Chern- 
Simons theory. The resulting energy functional is positive when evaluated for gauge fields 
which satisfy Gauss's law. Nair's expression for £ reads £ = £yu + £n, with 



v-D, G 



v-D, G T 



(4.19) 
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and involves the auxiliary matrix field G(x,\) G SU(N), with the property G(x,v) = 
Gt(z,-v). 

We show now that, in spite of formal differences, the energy densities in eqs. Q4.1& ) 
and ( |4. 19|) are, in fact, the same. To do that, we first replace D$ = (£)_ + D + )/2 and 
v • D = (£>_ - D+) /2 (where D+ = vD and D_ = v'-D). We thus get, instead of (|4~T9|) , 



8 g 2 



— Tr j 



D_, G 



£>_, G T 



£>+, G 



L> + , Gt]}. 



(4.20) 



By changing v — > — v in the second term inside the curly braces, we see that it gives the 
same contribution as the first one after angular integration. Furthermore, we note that 



the equations of motion for G(x, v), as written down in Ref. [|16[], are implicitly solved by 



a_ = G~ l A_G - ~G _1 a_G, 
9 

which relates G to our function a_ = v' • a. Writing this equation in the form 



(4.21) 



--G- 1 



£L, G 



G f 



G. 



(4.22) 



we conclude that eq. (|4.20| ) is equivalent to 



N 



3 2 
4 W ' 



47T 



Tr{( ( 



(4.23) 



where any direct reference to the field G has disappeared. It is then sufficient to recall that 
a_ — A_ = 2W° to see that as given by eq. (|4.23| ), is identical to our expression ( 4.18|) 
for Bind- Remark that in the arguments above we did not need the explicit expression of 
G. However, if one recalls that a_ = (i/(?)(<9_V A )V A_1 (see eq. ( |3.18j )), with V(x;v) defined 
by ( |2.13| ), one can easily verify that eq. ( |4.21| ) is solved by G(x, v) = V(x; v') V^(x; v). 



5 Conclusions 

The classical features of the HTL's reflect the long- wavelength, collective behaviour of the 
plasma particles, namely, their average motion over distances which are large compared to 
their mean separation. The coherent motion of the particles is described, to leading order 
in g, by density matrices obeying simple kinetic equations. In the abelian case, the kinetic 
equation is the linearized Vlasov equation and the elementary dynamics is essentially that 
of the Lorentz equation which governs the motion of a classical charged particle in an 
electromagnetic field. In the non-abelian case, the distribution functions are matrices 
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in color space and their time dependence involves in addition a "precession" in color 
space. We have shown that a remarkably simple gauge structure emerges from our kinetic 
equations once one assumes that the plasma particles are not significantly deviated by the 
background field. Remarkably, this gauge structure determines completely the collective 
dynamics of the particles at the order of the HTL's approximation. This explains why the 
requirement of gauge symmetry was, in fact, a sufficient condition to obtain the effective 
action generating the HTL's||], as well as the related field equations 0, p~3f| . 
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